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Abstract 



A metastable homogeneous state exists down to zero temperature in systems 
of repelling objects. Zero "fluctuation temperature" liquid state therefore 
serves as a (pseudo) "fixed point" controlling the properties of vortex liq- 
uid below and even around melting point. There exists Madelung constant 
for the liquid in the limit of zero temperature which is higher than that of 
the solid by an amount approximately equal to the latent heat of melting. 
This picture is supported by an exactly solvable large N Ginzburg - Landau 
model in magnetic field. Based on this understanding we apply Borel - Pade 
resummation technique to develop a theory of the vortex liquid in type II 
superconductors. Applicability of the effective lowest Landau level model is 
discussed and corrections due to higher levels is calculated. Combined with 
previous quantitative description of the vortex solid the melting line is located. 
Magnetization, entropy and specific heat jumps along it are calculated. The 
magnetization of liquid is larger than that of solid by 1.8% irrespective of 
the melting temperature. We compare the result with experiments on high 
T c cuprates Y ' Ba^Cu^O'j , DyBCO, low T c material (K, Ba)BiOs and with 
Monte Carlo simulations. 
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I. INTRODUCTION AND THE MAIN IDEA 



Abrikosov flux lines (vortices) created by magnetic field in type II superconductors 
strongly interact with each other creating highly correlated configurations like the vortex 
lattice. In high T c cuprates thermal fluctuations at relatively large temperatures are strong 
enough to melt the lattice. Several remarkable experiments demonstrated that the vortex 
lattice melting in high T c superconductors is first order with magnetization jumps [1] and 
spikes in specific heat [2]. Magnetization and entropy jumps were measured using local 
Hall probes [1], SQUID [3,4], torque magnetometry [5,6] and integrating the specific heat 
spike [2,7]. It was found that in addition to the spike there is also a jump in specific heat 
in YBCO which was measured as well [2,7,8]. These precise measurements pose a question 
of accurate quantitative theoretical description of thermal fluctuations in vortex matter. 
The melting line in high T c cuprates has been studied mainly not very far from T c . In this 
part of the phase diagram the Ginzburg - Landau (GL) approach is generally appropriate 
to describe thermal fluctuations near T c [9,10]. The GL model is however highly nontrivial 
even within the lowest Landau level (LLL) approximation valid at relatively high fields. 
This simplified model has only one parameter: the dimensionless LLL scaled temperature 

(T - T m f(H))/(TH) 2 / 3 (defined precisely in eq.(13a) below). Over the last twenty 
years great variety of theoretical methods were applied to study this model. Brezin, Nelson 
and Thiaville [11] applied the renormalization group (RG) method on the one loop level 
description of the vortex liquid. No nontrivial fixed points of the (functional) RG 
equations were found and they concluded therefore that the transition from liquid to the 
solid is first order [12]. Another often used approach applicable also beyond the range of 
validity of GL model is to use elasticity theory description of the vortex lattice and 
Lindermann criterion to determine the location of melting line [13]. However all those 
approaches do not provide a quantitative theory of melting since these are one phase 
theories and in order, for example, to calculate discontinuities at first order transition 
accurate description of both phases is necessary. 

Two perturbative approaches were developed and greatly improved recently to describe 
both the solid and the liquid phases in the LLL GL model. The perturbative approach on 
the liquid side was pioneered long ago by Thouless and Ruggeri [14]. They developed a 
perturbative expansion around a homogeneous (liquid) state in which all the "bubble" 
diagrams are resumed. Unfortunately they found that the series are asymptotic and 
although first few terms provide accurate results at very high temperatures, the series 
become inapplicable for less than —2, which is quite far above the melting line (believed 
to be located around ~ —10). We recently obtained the optimized Gaussian series [15] 
which are convergent rather than asymptotic with radius of convergence of a*r = —5 still 
unfortunately above the melting point. 

On the solid side, long time ago Eilenberger and Maki and Takayama [16] calculated the 
fluctuations spectrum around Abrikosov's mean field solution. They noticed that the 
vortex lattice phonon modes are softer than that of the acoustic phonons in atomic 
crystals and this leads to infrared (IR) divergences in certain quantities. This was initially 
interpreted as " destruction of the vortex solid by thermal fluctuations" and the 
perturbation theory was abandoned. However the divergences look suspiciously similar to 
"spurious" IR divergences in the critical phenomena theory and recently it was shown that 
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all these IR divergences cancel in physical quantities [17]. The series therefore are reliable, 
and were extended to two loops, so that the LLL GL theory on the solid side is now 
precise enough even around melting point. 

Therefore the missing part is a theory in the region — 10 < a T < — 5 on the liquid side. 
Moreover this theory should be very precise since free energies of solid and liquid happen 
to differ only by few percents around melting. Closely related to melting is the problem of 
nature of the metastable phases of the theory. While it is clear that the overheated solid 
becomes unstable at some finite temperature, it not generally clear whether the overcooled 
liquid becomes unstable at some finite temperature (like water and other molecular liquids, 
which however has a crucial attractive component of the intermolecular force) or exists all 
the way down to T = as a metastable state. The Gaussian (Hartree - Fock) variational 
calculation, although perhaps of a limited precision, is usually a very good guide as far as 
qualitative features of the phase diagram are concerned. Such a calculation in the liquid 
was performed long ago [14], while a significantly more complicated one sampling also 
inhomogeneous states (vortex lattice) was obtained recently [18,19]. The gaussian results 
are as follows. The free energy of the solid state is lower than that of the liquid for all 
temperatures lower than melting temperature a™. The solid state is therefore the stable 
one below a™, becomes metastable at somewhat higher temperatures and is destabilized at 
ax = —5. The liquid state becomes metastable below the melting temperature, but unlike 
the solid, does not loose metastability all the way down to a™ = — oo (T = 0). The 
excitation energy of the supercooled liquid approaches zero as a power e ~ l/a T . This 
general picture is supported in section III by an exactly solvable large N Ginzburg - 
Landau model of vortex matter in type II superconductors. 

Meantime in different area of physics similar qualitative results were obtained. It was 
shown by variety of analytical and numerical methods that liquid (gas) phase of the 
classical one component Coulomb plasma exists as a metastable state down to very low 
temperature possibly T = [20] . The quantum one component plasma - electron gas also 
shows similar feature [21]. One considers quantum fluctuations instead of the thermal ones 
(an analog of inverse temperature is coupling r s , see section III). It seems plausible to 
speculate that the same phenomenon would happen in any system of point like or line like 
objects interacting via purely repulsive forces. In fact the vortices in the London 
approximation are a sort of repelling lines with the force even more long range than 
Coulombic. This was an additional strong motivation to consider the above scenario in 
vortex matter. In section III we provide both theoretical and phenomenological evidence 
that the above scenario is the correct one. 

Assuming absence of singularities on the liquid branch allows to develop an essentially 
precise theory of the LLL GL model in vortex liquid (even including supercooled liquid) 
using methods of theory of critical phenomena [22,23]. The generally effective 
mathematical tool to approach a nontrivial fixed point (in our case at zero temperature) is 
the Borel - Pade (BP) transformation [23]. Before embarking on this program in following 
sections, we address several subtleties which prevented the use and acceptance of the BP 
method in the past and then combine it with the recently developed LLL theory of solids 
to calculate the melting line and the magnetization and the specific heat jumps across the 
line. Very early on Ruggeri and Thouless [14] tried to use BP (unfortunately a 
"constrained" one, so that it interpolates smoothly with the solid, assumption we believe is 
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incorrect) to calculate the specific heat without much success. It was shown by Wilkin and 
Moore [24] that the constrained BP does not converge, while the results for unconstrained 
BP were inconclusive. They attributed this to the limited order of expansion known at 
that time. Subsequent attempts to use BP for the calculation of the melting line using 
longer series also ran into problems. Hikami, Fujita and Larkin [25] tried to find the 
melting point by comparing the BP energy with the one loop solid energy and obtained 
(It = —7. However their one loop solid energy was incorrect and, in any case, it was not 
precise enough (as will become clear below the two loop contribution cannot be neglected). 
The LLL GL model was also studied numerically in both Lawrence - Doniah model (a 
good approximation of the 3D GL for large number of layers) [26,27] and in 2D [28] and by 
a variety of nonperturbative analytical methods. Among them the density functional [29] , 
1/N [30-32], dislocation theory of melting [33] and others [34]. 

As we show in this paper, the BP liquid free energy combined with the correct two loop 
solid energy computed recently gives scaled melting temperature a™ = —9.5 and in 
addition predicts other characteristics of the model. The magnetization of liquid is larger 
than that of solid by 1.8% irrespective of the melting temperature (the specific heat jump 
is about 6% and decreases slowly with temperature in YBCO). Brief account of these 
results was published [18]. 

In addition to the theory of melting, we considered overcooled liquid and calculated 
magnetization and specific heat curves. Since the metastable overcooled liquid state exists 
all the way down to zero temperature, we can define liquid Madelung energy. Looking at 
the melting process from the low temperature side for both the liquid and the solid we find 
that the Madelung energy of liquid is larger than that of the solid approximately by the 
latent heat of melting. Our magnetization curves agree quite well with Monte Carlo 
simulations of the LLL GL [26] , and almost perfectly for specific heat in 2D by Kato and 
Nagaosa in ref. [28]. 

In addition we consider in this paper several "phenomeno logical" issues, some matter of 
significant disagreement. First is the range of applicability of the LLL model. We find that 
in order to describe experimental reversible magnetization of YBCO at lower fields, higher 
Landau levels (HLL) corrections should be incorporated. We therefore clarify in section V 
the role of the HLL modes. Experimentally it was claimed that one can establish the LLL 
scaling for fields above 3T [35] . A glance at the data however shows that in normal state 
(above T c ) the LLL scaling for magnetization curves is generally very bad. Most of the 
HLL effects can be taken into account by just renormalizing parameters of the LLL model. 
Therefore one should use the "effective LLL" in which HLL were "integrated out". To 
clarify this often salient feature we explicitly perform this integration within a self 
consistent approach in section V. It was noted by Koshelev [36] and others that, to 
calculate magnetization, one has to carefully account for renormalization of the free energy 
since it is field dependent. Then we calculated the leading correction the effective LLL and 
compared with experiments. It is found that although the LLL contribution to 
magnetization is much larger that the experimentally observed one above T c , it is nearly 
cancelled by the HLL contributions. This explains the breaking of the LLL scaling in the 
normal state. 

The paper is organized as follows. The model is defined and its applicability range 
discussed in section II. In section III, supercooled liquid in general repulsive system will be 
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discussed. We provide evidence for the scenarios outlined above using both a solvable (for 
a large number of components) Ginzburg - Landau model and analyzing numerical results 
in other systems. In section IV the LLL model is solved and the melting theory of vortex 
lattice is presented and compared to experiments. In section V, the HLL corrections are 
discussed and the magnetization curves are compared with experiments. 
Phenomenological issues are addressed in sections IIB (assumptions), IVC (melting line, 
Ginzburg parameter fit for various materials), IVD (magnetization, entropy jumps), IVE 
(specific heat jumps) and VD (reversible magnetization curve), so readers not interested in 
theoretical details can directly proceed to these sections. 

II. THE GL MODEL AND ITS BASIC ASSUMPTIONS. 

A. The GL model 

On the mesoscopic scale a 3D superconducting materials with not very strong asymmetry 
along the z axis are effectively described by the following Ginzburg-Landau free energy 
functional: 

F[W, A] = J d*x^- W\ 2 + ^\d-M 2 - a(TM\ 2 + b -m 4 + (1) 

involving the order parameter field ip and magnetic field B. The external constant 
magnetic field is described by the vector potential in Landau gauge A = (Hy, 0, 0). The 
covariant derivative is defined by D = V — 2niA/® o ,& = hc/e*(e* = 2e). The 
microscopic thermal fluctuations are integrated out and, as a consequence, coefficients a, b' 
and m depend on temperature. Mesoscopic thermal fluctuations of the order parameter 
are described by the partition function: 

^/xWPAexpj-^f^I} (2) 

Our aim is to quantitatively describe the effects of thermal fluctuations of high T c cuprates 
of the YBCO type and some other "strongly fluctuating" 3D materials. 

B. Assumptions 

The use of the above GL energy hinges upon several physical assumptions. They are listed 
below. 

(1) Continuum 3D model 

We use the anisotropic GL model despite the well established layered structure of the high 
T c cuprates for which models of the Lawrence - Doniah type are more appropriate. Effects 
of layered structure are dominant in BSC CO or Tl compounds (anisotropy very large: 
7 = yj m c / m ab > 1000) and noticeable for cuprates with anisotropy of order 7 = 50 like 
LaBaCuO, strongly underdoped YBCO (see however [37]) or Hgl223. The requirement, 
that the 3D GL can be effectively used, therefore limits us to optimally doped YBCOj^s 
(or slightly overdoped or underdoped) for which the anisotropy parameter is not very large 
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7 = 4 — 8 [38]), DyBCO and possibly Hg\T2\ which has a slightly larger anisotropy. 
There is no such problem in recently discovered isotropic "fluctuating" superconductor 
(K, Ba)BiO ?> [39]. 

(2) Range of validity of the mesoscopic (GL) approach 

The GL approach generally is an effective mesoscopic approach. It is applicable when one 
can neglect higher order terms in the functional eq.(l) typically generated, when one 
"integrates out" microscopic degrees of freedom. The leading higher dimensional terms we 
neglect (as "irrelevant") are and higher (four) derivative terms like |D 2 ?/>| 2 . This 
naively leads to a condition that 1 — t — b is smaller than 1. Here and in what follows 

t = T/T c ; b ee H/H c2 » B/H c2 . (3) 

The applicability line 1 — t — b < 0.2 for YBCO is plotted in Fig. 1. We also will consider 
a model invariant under rotations in the ab plane. Noninvariant models sometimes can be 
rescaled it to m a ~ m& = m a 5 [10]. For several physical questions those assumptions are not 
valid because neglected "irrelevant" terms might become "dangerous". For example the 
question of the structural phase transition into the square lattice is clearly of this type 
[40]. It is known that even assuming m a /m b = 1 in low temperature vortex lattices in 
YBCO, rotational symmetry is broken down to the fourfold symmetry by the four 
derivative terms. However there is no significant correction to, for example, the 
magnetization from those higher dimension terms. 

(3) Expansion of parameters around T c 

Generally parameters of the GL model of eq.(l) are complicated functions of temperature 
which are determined by the details of the microscopic theory. We expand the coefficient 
a{T) near T c : 

a(T) = T c [a(l -t)- a {I - t) 2 + ...]. (4) 

The second and higher terms in the expansion are omitted and therefore, when 
temperature deviates significantly from T c , one cannot expect the model to have a good 
precision. We note that recently measured H C 2(T) is linear in T in a wide region near T c in 
both YBCO and (K,Ba)Bi0 3 [41,39]. 

(4) Constant nonfiuctuating magnetic field 

For strongly type II superconductors like the high T c cuprates not very far from H C 2(T) 
(this easily covers the range of interest in this paper, for the detailed discussion of the 
range of applicability beyond it see ref. [42]) magnetic field is homogeneous to a high 
degree due to superposition from many vortices. Inhomogeneity is of order 1/k 2 ~ 10 -3 . 
Since the main subject of this study is thermal fluctuation effects of the order parameter 
field, one might ask whether thermal fluctuations of the electromagnetic field should be 
also taken into account. Halperin, Lubensky and Ma considered this question long time 
ago [43]. The conclusion was that they are completely negligible for very large k. Upon 
discovery of the high T c cuprates, the issue was reconsidered [44] and the same result was 
obtained to a very high precision. Therefore here magnetic field is treated both as constant 
and nonfiuctuating (B = H) and the last term in eq.(l) can be omitted (to precision of 
order l//t 2 ). However when we calculate the magnetization, M = (B — H)/4tt which is of 
order 1/k 2 , higher order correction must be considered. 
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Recently it was claimed that the "vortex loops" fluctuations are important and even might 
lead to additional phase transition at field of order Gi H c2 [45]. This is of order 100G for 
the materials of interest listed in Table 2 and therefore is irrelevant for physics discussed in 
this paper. Note that Gi in the papers discussing the vortex loops [46] physics is assumed 
to be much larger. We discuss this issue in section IVC 
(5) Disorder. 

Point - like disorder is always present in YBCO. For example magnetization becomes 
irreversible. Melting line of the optimally doped or underdoped samples bends towards 
lower fields [7] and signs of the second order transition appear at 12T [47]. However in 
some samples like fully oxidized YBa 2 Cu 3 7 [6] and DyBa 2 Cu 3 7 [8,48] the disorder 
effects are minor especially at temperature close to T c . In the maximally oxidized YBCO 
[6] the second order transition associated with disorder is not seen even at highest available 
fields (30T). Certain aspects of the disorder problem were addressed in the framework of 
GL theory [49], elasticity theory [50] and phenomenological approach based on the 
Lindermann criterion [13]. 



C. Landau level modes in the quasi momentum basis 

Assuming that all the requirements are met we now divide the fluctuations into the LLL 
and HLL modes. Throughout most of the paper will use the coherence length 
£ = \Jh 2 / (2m ab aT c ) as a unit of length, T c as unit of temperature and dHc ^ Tc ^ T c = ^fp- as 
a unit of magnetic field. As we mentioned above, we assume constant magnetic induction 
B = hH c2 which is slightly different from the external magnetic field H = hH c2 . After 
rescaling eq.(l)by x — > £x, y — > £y, z — > ip 2 — > ^f^ip 2 (7 = i/m^/nu) one obtains the 
Boltzmann factor: 



J T 00 



<TX 



l -W\ 2 + \\dM 2 



a h + - 



| 2 + ^| 4 + ^ 



h) 5 



(5) 



where dimensionless parameter 



\p2GmH 



(6) 



characterizes the strength of thermal fluctuations. The commonly used dimensionless 
Ginzburg number is defined by 



Gi = 



1 /327re 2 /t 2 £T c 7 

2 V ~&h 2 



(7) 



And 



an = 



(8) 



defines the distance from the mean field transition line. It is convenient to expand the 
order parameter field in a complete basis of noninter acting theory: the Landau levels. In 
the hexagonal lattice phase the most convenient basis is the quasi- momentum basis: 
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if>(x,y,z) 



%/2(2tt) 



27r ) ^ n=0 



(9) 



Here is the eigenstate of the n th Landau level e n = (n + 1/2)6 with two dimensional 

quasi- momentum k with the hexagonal symmetry: 



2"- 1 n!a A ^ "'^ ' " ' Vb a A 

i=— 00 



x exp < « 



7T 



'(»-!) , i^^-jt) , Tfc 



(10) 



where a A = y The function ^ = y?k=o describes the Abrikosov lattice solution [9]. 

Even in the liquid state which is more symmetric than the hexagonal lattice, we find it 
convenient to use this basis. 

Naively, if the magnetic field is sufficiently high, the energy gap of the order b separating 
the n = LLL modes from the HLL is very large it is reasonable to keep only the LLL 
modes in eq.(5). The dominance of the LLL modes for melting was discussed in ref.( [11]), 
and Pierson and Vails in ref. [35] and we will discuss it in more detail in section V. In the 
rest of this section, we consider the LLL GL model. 



D. The LLL scaling 

Using the LLL condition |D?/>| 2 = fe^l 2 , the free energy simplifies: 



/ = - / d 6 x 



\\dM 2 -a h \^\ 2 + \\^+ K2{h A ^ 



(11) 



There is no longer a gradient term in directions perpendicular to the field and 
consequently the model possesses the LLL scaling [51]. After additional rescaling 



x — > x/yb, y — > y/vb, z — > z ( j , ip — > ( -^j^ J 4>, the dimensionless free energy 
takes a form: 



/ 



Airy/2 



(Fx 



( boo \ ~ 4/3 (b - hf 



(12) 



Minimizing it with respect to b leads to magnetization b — h of the order ■ This means 
that in the strongly type II limit (k » 1) the last term is of the order 1/k 2 and can be 
neglected. The theory has a single dimensionless parameter, the Thouless scaled 
temperature defined by: 



bui 



2 5/2 



7T 



-2/3 



a h . 



(13a) 
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The Gibbs free energy density in the newly scaled model is defined as: 

9M = log J DijDij* exp {-/M} ; (14) 

which is also a function of ar only (47r\/2 is the scaled "temperature".). The relation to 
the original Gibbs free energy is: 



III. OVERCOOLED LIQUID AND THE T=0 FIXED POINT OF THE LLL 

MODEL 

A. Mean field approximation: absence of the finite temperature pseudo critical point 

for the vortex liquid. 

The energy of the hexagonal solid in mean field (neglecting mesoscopic thermal 
fluctuations) is [9]: 

n 2 H 2 
sol _ a T . r sol _ n c2 2 / 1fi \ 

where Pa = 1-1596 and the subscript "M" underlies similarity to the Madelung energy of 
atomic solids. The major fluctuations contribution to the solid free energy is due to the 
"phonon" modes. In harmonic approximation it is proportional to the fluctuation 
temperature T = a T 3 ^ 2 : 

gf = 2.848 \a T \ 1/2 ; Gf = CfT; (17) 
C f = 2.848^ 

At low fluctuation temperatures one can neglect the T dependence of at — —(1 — b)/2. 
Solid becomes unstable at ar = —5 according to the self consistent (gaussian) 
approximation [19]. 

In the (homogeneous) liquid state order parameter vanishes and the contributions to free 
energy come solely from fluctuations. The gaussian ("mean field") approximation to the 
free energy [14] is 

4 = 4^-4/5, (18) 
where the excitation energy e is given by a solution of the cubic "gap equation" 

e 3/2 - a T Ve ~ 4 = 0. (19) 

The liquid state becomes metastable below the melting temperature, but unlike the solid 
above melting, does not loose metastability at certain "spinodal" point [52]. It persists all 
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the way down to T = 0. The excitation energy of the supercooled liquid approaches zero as 
a power e ~ 16/ a 2 -. For a? — ► — oo, the scaled energy eq.(18) has an expansion in 
l/a\ oc T 2 for small fluctuation temperature T (the radius of convergency of the expansion 
extending to a T = —3). Therefore the liquid despite having energy larger than that of solid 
becomes (pseudo) critical [53] at zero temperature. Physical quantities "around" this point 
exhibits a power behavior with characteristic (pseudo) critical exponents. The metastable 
liquid state has a distinct Madelung energy 

<' = -«§4. (20) 

As temperature increases the difference between the solid and the liquid becomes smaller 
and vanishes at melting. Generally one expects a linear correction at small T: 

G Uq = G l lf + C[ iq T. (21) 

Since the expansion of the mean field free energy is in T 2 : C\ q = 0. Comparing the solid 
free energy eqs. (16,17) with eq.(21), we get the melting temperature a™ = —6.3. We 
therefore conclude that in this approximation the supercooled liquid state exists down to 
its pseudo critical point at zero temperature. Moreover the pseudo critical point might 
govern the behavior of the liquid phase to temperature as high as the melting point. 



B. The large N approximation and the Lopatin - Kotliar model of the Abrikosov 

lattice melting 

It is important to confirm the above scenario in an exactly solvable model. The simplest 
model of this kind is the multicomponent GL model. The LLL GL theory can be 
generalized (in several different ways) to an N component order parameter field 
r,a = l,...,N: 



f = 7= / d A x 



(22) 



The large N limit of this theory can be solved in a way similar to that in the N component 
scalar models widely used in theory of critical phenomena [22]. The simplest case v — 1 
has been considered in ref. [30]. It was found that the homogeneous state is stable at all 
temperatures. Under assumption that the conventional Abrikosov lattice takes over at low 
temperatures it supported the original conjecture by Brezin et al. [11] that melting of the 
flux lattice is a first order phase transition. However it was shown (by explicit numerical 
evaluation) in [12] that the low temperature ground state in that model is not the 
Abrikosov lattice state in which just one component has a nonzero expectation value 
(similar to the one component Abrikosov lattice). The "true" ground state has infinite 
degeneracy. Different ground states at large N are markedly different from the hexagonal 
lattice. The case v = 2,. in which the Abrikosov lattice state is a stable ground state, was 
first introduced in [32] (in what follows it will be referred to as the LK model). Eq.(22) is a 
slight generalization including both models studied in ref. [30,32]. We find that in fact all 
models with v > 2 possess such a stable lattice state. 
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A straightforward method to develop the 1/N expansion with the last component of ip N 
having the expectation value oc ipA = y?k=0' describing the hexagonal lattice (see see 
eq.(10)), is to shift this field tp N (x,y,z) — > tp N (x,y,z) + y/Nc(pA(x,y), where c is a (real) 
constant. Then one introduces Hubbard - Stratonovich (HS) fields p, x [32] via free energy: 

nr, P , x\ = ^= (\\w a \ 2 + iyp + arW\ 2 + *«?\<PA\ 2 m 2 + ^ [{<?<& + xw b r b + «)] 

N I v 1 — v . l9 \ „ . . 

^ -P 2 H ^Ixl 2 > +JV /„/ + ... (23) 



Atx^2\2' 2 
Here the " nonfluctuating part" is the Abrikosov free energy density 



ff = — 



2 , Pa 4 

a^c H c 

2 



(24) 



We omitted several cubic terms which do not influence the leading order in 1/N. 
Integrating over the fluctuating the ip a fields one obtains the effective scaled Gibbs energy 
density (the calculation is very similar to that in [19], where technical details can be found): 

if = arc? + f c> - + ^Ixf ) i + 2 {rftfT) + ■ (25) 

The spectrum has two branches: 

eo,A(k) = a T + v ((3 k c 2 + p k ) ± | (1 - v) (c 2 lk + Xk)\)- (26) 

To have a stable perturbative Abrikosov solution, the spectrum should be positive definite 
for Pk = Xk = 0. Thus we demand —v/2 + {y — 1) > or v > 2, as stated above. Here 
both HS fields 

Pk = (p(x)\(p k (x)\ 2 ) x ; Xk = (x(x)(p* k (x)(fi*_ k (x)) x (27) 
and the constant c should minimize free energy g e ff- 

C. The inhomogeneous (solid) solution 

The minimization with respect to p (x) and x ( x ) leads to 

P(x) = (|^(^)| 2 {[ e o(k)]- 1/2 + [^(k)]- 1/2 }) fc (28) 

sign (I - v)x(x) = ( fk(x)ip- k (x) ° lk ~j~ Xk {[e (k)] _1/2 - {e A (k)]~ 1/2 j 
\ \c ik + Xk\ >- J 

which, in terms of Fourier harmonics of the hexagonal lattice, takes a form: 

Pi = (Pi-k {[eo(k)]- 1/2 + [e A (k)]- 1 / 2 }) k (29) 
signal -v)n = U, ^ 1 {[^(k)]- 1/2 - Mk)]- 1/2 } 
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The lattice functions Pk,1k,1k,i are defined in appendix of [19]. The only consistent 
solution preserving the hexagonal symmetry is Xk = Xdk, and the above equation will 
simplify to: 

sign (1 - u) XcPa = (sign (c 2 + Xc ) % { [e (k)]- 1/2 - [e A (k)]' 1/2 }) fc . (30) 

For the LK model [32], v = 2,this leads to: Xc (3 A = { Mk)]~ 1/2 - [e (k)]" 1/2 |^ . 
Finally the set of the minimization equations ( X > 0) is 

= a T + (3 A c 2 + 2 (ft { [e G (k)]- 1/2 + Mk)]- 1/2 }) ^ + 
(, fc {[ eo (k)]- 1 / 2 -[ e ,(k)]- 1 / 2 }) k 
XcPa = (vk {[e A m^ 2 - [ eo (k)]- 1/2 }) k (31) 



Pl = (^k{[eo^)r 1/2 + [e A (k)r 1/2 } 



and 



e A k ) = ar + 2f3 k c 2 + 2p k ± (c 2 + Xc ) 7fc . (32) 
Following formulas will be useful for the calculation of the free energy: 

(p 2 ) = (a_* {Mk)r l/2 + Mk)]-/ 2 } {Mir 172 + Mi)r 1/2 }) M 

<|X| 2 > =^[(^{^(k)]- 1 / 2 -[ eo (k)]- 1 / 2 })J 2 (33) 

The equations in eq.(31) can be solved using mode expansion [32,19]. The spectrum can be 
written as follows 

e (k) = E(k) + A Vk , e A (k) = E(k) - A Vk , (34) 
with E(k) expanded in modes 

E(k) = Y,E n (3 n (k), (35) 

where 

oo 

(3 k = J2 e M-^n/V3](3 n (k), f3 n (k)= ex P[* k - x ]- (36) 

The integer n determines the distance of a points on reciprocal lattice from the origin. The 
effective "expansion parameter" is exp[— 27r/v^3] = 0.0265 and coefficients decrease 
exponentially with n [19] as can be seen from Table 1. 
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Table 1 

Coefficients of the mode expansion for the solid solution 



ax 


9 


Ex 


Ei 


E 3 


A 


-4.6179 


-3.43164 


0.728715 


-0.0022412 


-0.00001227 


0.6167 


-5 


-4.96636 


1.92669 


0.0717767 


0.00003881 


2.0331 


-10 


-34.3165 


6.29543 


0.355908 


0.00023872 


7.2718 


-20 


-159.826 


13.8477 


0.842385 


0.00058357 


16.3036 



The solution disappears at ar = —4.6179. At this point the solid is no longer a metastable 
state, eyi(k) is a gapless mode and e^k) — > const.k 2 for k — > 0. For perturbative 
spectrum, e^(k) — > const.ls 4 for k — > 0. 



D. Melting in the LK model 

The energy corresponding to the solid solution of the minimization equation eq.(31) 
calculated from 

9 if = a ^ 2 + Y c4 -( p2 - 1 2 lxl2 ) +2(v / ^(kT+V / ^(k)) k (37) 

is given in Table 1. The convergence of the mode expansion is exponential. 

For liquid, we impose rotationally invariant Ansatz with c 2 = 0, x — and obtain the gap 

equation 

, = (38) 

which minimizes energy 

9u g = -P 2 + W a r + 2p- (39) 

The results for both the liquid and the solid free energy are ploted on Fig. 2. The melting 
point appears at ax = —5.15. 

It is well approximated in the whole region by its low temperature expansion in powers of 
l a T|~ 3 ^ 2 (which is proportional to the "fluctuation temperature" T assuming that at low 
temperatures ~ — (1 — 6)/2) 

l— = c ^ + c fT + cfT 2 ...,T = |a T |- 3/2 , (40) 



CLrp 



c s ° l = — ; cf = 2.84835; cf = -2.54087. 

2pa 

The first two terms are the same as for the one component model, while the two loop 
correction is different. 

Similarly the liquid energy can be expanded, but this time in powers of square of the 
"fluctuation temperature T 
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il = c g + c f T + 4 q T 2 ...,T = |a T |" 3 / 2 , (41) 



liq _ . /i<j _ n . Jig _ fi . Liq _ 9 n 
C M — 7> C l,3,.. — U J C 2 _ °i c 4 _ zu - 



Here the first term is the "Madelung energy" of liquid at zero fluctuation temperature. 
Note that, as in the mean field approximation to the one component theory, there is no 
term linear in T (the harmonic approximation). This means that the specific heat vanishes 
at zero temperature. Retaining just the Madelung and the harmonic term for solid we 
estimate the melting temperature in the linear approximation: 

sol _ Hq 

T = Cm Cm 

r hq - r s ° l 
C\ c 1 

The latent heat heat in the same approximation is: 

AC/ = dg - eg. (43) 

Numerically this melting temperature T m = 0.064 corresponding to ax = —6.25 and the 
latent heat AU = 0.18 should be compared with the exact results: T m = 0.086 
(or = -5.15), AU = 0.122945. 

To conclude obtains the first order melting. Supercooled liquid persists as a metastable 
state all the way to zero temperature. We emphasize that this means that the matching of 
the (Borel - Pade approximant to) liquid to solid energy at T = employed in [14] to 
improve convergence of the series is not only ineffective [24] , but should lead to an incorrect 
result. Liquid and solid energies are different in the limit of zero fluctuation temperature. 



E. General hypothesis about melting of lattices made of repelling objects 

In atomic liquids, an attractive long range force is generally present. As a result the 
supercooled liquid state loses its metastability at an end point (spinodal) [52]. Lovett 
argued on general grounds long time ago [54] (stability analysis of approximate set of 
relations between density correlators) that for certain purely repelling interactions the 
spinodal point disappears (shifted to zero temperature) and is recovered when the 
attractive interaction is introduced. The existence of a metastable overcooled liquid down 
to zero temperature for repelling particles therefore might be quite general. The best 
studied example of mutually repelling particles is the classical one component Coulomb 
plasma. We assume that, as in the vortex system at low temperature, the supercooled 
liquid has a Madelung energy and moreover its free energy has a low temperature 
expansion. The free and internal energies at low temperatures can be expanded as: 

F sol,liq = C sol,liq + C sol,liq Tf + ^liqrpj _ 3^ log T/ + (44) 
Q 

TjsolMq ,rysol,liq . u rp s^sol,liqrp2 , 

u = c M +-J/-c 2 lf + ..., 

where scaled temperature Tf is inverse of to the dimensionless plasma parameter 

r = (^p-) y f° r density n s . Existent Monte Carlo simulations of the internal energy in 
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the stable and metastable region of the 3D one component Coulomb plasma [55] can be 
well fitted (see solid line on Fig. 3) by 

C 1 ^ = -0.89186, C l 2 iq = -23.89. (45) 

It underestimates the internal energy at higher temperatures. Note that this fit is quite 
different from a variety of the fractional power expressions used at higher temperatures. 
One of the more successful (not very far from the melting point) liquid theories based on 
density functional approach is [56]: 

U Uq = -0.9 + 0.388 Tf 5 (46) 

(dashed line in Fig. 3). Note however that these expression cannot be continued to T = 0, 
since finite packing parameter is assumed. Unfortunately the coefficient C\ 9 cannot be 
deduced from internal energy only, while free energy is not available at large coupling. 
In the solid the analytical calculation gives [57] 

C$ = -.895929, Cf = -1.8856, Cf = -10.84. (47) 

This was corroborated by recent simulations [58]. Using the linear in T approximation 
eq.(43) we obtain latent heat AU = 0.0041, which should be compared with MC 
simulation result [59] AU = 0.043. From the measured melting temperature [57] 
T rn = 1/172 and coefficients fitted above one deduces 

C l l q = -2.5. (48) 

One observes that even at melting the linear approximation is justified (the C^ 9 
contribution to free energy account for less than 10% of the linear one). It would be very 
interesting to simulate the 3D Coulomb plasma at even stronger coupling T > 200 to verify 
the existence of expansion of supercooled liquid free energy as in eq.(44).We expect that 
other liquids with purely repulsive interactions like the Yukawa (screened Coulomb) 
studied recently in connection with "dusty plasma" or colloid suspensions physics [60] or 
even hard core repulsion lead to qualitatively similar result. An intriguing question is 
whether structure function is universal in the zero temperature limit of the liqiud phase. 
Since ideal liquid is (pseudo) critical, certain universal properties are expected. 
Even more closely related to the vortex system is the quantum one component plasma - 
electron gas. The quantum particle is described by a fluctuating line very analogous to a 
thermally fluctuating vortex line. Here we consider 2D electron gas at zero temperature 

ff = -SS^« + 2^1^77- (49) 



The path integral of this system is quite analogous to a system of repelling flux lines. The 
quantum fluctuations in 2DEG replace the thermal fluctuations of the vortex system. In 
quantum partition function 



D^* exp 



(50) 
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with the fermionic (Grassmannian) field \I/ replacing the bosonic field in eq. (2). The 
(Wigner) solid solution can be well approximated by the expansion in quantum 

fluctuations T q = l/^T s = h^^- [21], 

where C$ = -2.2122, Cf = 1.6284, Cf = 0.058. For unpolarized liquid, a very good fit 
in wide range of densities is [21]: 

E u q = _ 8^2 a (T q + oQ 

3tt 9 T| + fll T2 + a 2 T q + a 3 
a = -0.3568, a x = 1.13, a 2 = 0.9052, a 3 = 0.4165. 

The fit in the low tempeature region 

E liq = + C[ ig T q + ... 

gives liquid Madelung energy C l $ = -2.18154, while leading correction is C[ iq = 1.45266. 

The transition to Wigner crystal occurs at T™ = 0.174 corresponding to 

r s = 33. Variational MC simulation [21] indicates that the transition occurs at r s = 37. 



IV. BOREL - PADE METHOD APPLIED TO THE LLL MODEL. MELTING 
LINE, MAGNETIZATION AND SPECIFIC HEAT. 

A. The BP method applied to liquid energy 

As we have seen above, within mean field the liquid branch exhibits a pseudo critical point 
[53] at T = 0. It is well known that in the theory of critical phenomena one can obtain an 
accurate description in the critical region by applying the Borel - Pade methods to 
perturbation expansion at "weak coupling" [23]. In technical terms there exists a 
renormalization group flow from the weak coupling fixed point towards the strongly couple 
one [22]. We therefore start with (the renormalized) weak coupling (high temperature or 
non-ideal gas) expansion. 

The liquid LLL (scaled) free energy is written as [14] 

g aq = 4e 1 ^ 2 [l + h(x)]. (51) 
The function h can be expanded as 

h(x)=Y,CnX n , (52) 

where the "small parameter" x = |^~ 3//2 is defined as a solution of the Gaussian gap 
equation for the excitation energy e,eq.(72). The coefficients c n can be found in [25]. The 
consecutive approximants are plotted on Fig. 4 as dashed lines (Tl to T9, TO being 
equivalent to the Gaussian mean field). One clearly sees that the series are asymptotic and 
can be used only at > —2. One can improve on this by optimizing the variational 
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parameter e at each order instead of fixing it at the first order calculation. The procedure 
is rather involved, see [61], however the optimized perturbation series is convergent with 
radius of convergence about = —5 (see dash dotted lines 1 to 9 on Fig. 4). Now we 
construct the BP series and compare them with the optimized perturbation series results 
for ax > — 5 . 

We will denote by hk (x) the [k, k — 1] BP transform [23] of h(x) (other BP approximants 
violate the correct low temperature asymptotics). The BP transform is defined as 

roc 

hk — I h k (xt) exp (— t) dt (53) 
Jo 

where hk is the [k, k — 1] Pade transform of Y^n=i £n T~> namely a rational function 

y^fe a . x i 

^fc-i 6 ' i with the same expansion at small x as the original function. 

For k = 4 and k — 5, the liquid energy converges to required precision (0.1%), see Fig. 4. 
On this Figure only k — 3 and 5 are shown since k = 4 practically coincides with the 
latter. In what follows we will use h 5 as the best available approximation of the liquid 
branch. The liquid energy completely agrees with the optimized Gaussian expansion 
results [15] till its radius of convergence at ax = —5. We therefore conclude that k = 5 is 
quite sufficient for our purposes. 

Since the metastable liquid state exists at all temperatures one can consider the T = 
limit. One finds: 

q lig (ar) , s 

- [ — ► 0.964 54 

g sol {a T ) v ; 



.,2 



for ax — > — oo. For g sol (aT), the leading term in this limit is — which is the Madelung 

2 

energy of the solid. The leading term for g hg (ar) is —0.964^-. Usually the Madelung 
energy for the solid phase of the point particle system is realized by minimizing the 
potential energy of the system (the minimum is often obtained by taking the hexagonal 
lattice for the repulsive system in 2D). In this vortex system, we can have the supercooled 
liquid down to ar — > — oo (T — > 0). The leading term for the overcooled liquid energy, or 

a 2 

the Madelung energy of the liquid is therefore equal to —0.964^-, which is slightly larger 
than the Madelung energy of the solid. This limit, the "ideal liquid", however cannot be 
thought as a minimization of a potential energy. 



B. Melting line. Comparison with Monte Carlo simulations and Lindemann criterion. 

The solid energy to two loops is [17,19]: 

9s* = - J 1 + 2-848 \a T \ 1/2 + — . (55) 
2j3 A a T 

where (3a = 1.1596. On Fig. 1 of ref. [18] we plot the energies of solid and liquid. They are 
very close near melting (see the difference on inset of this figure). We find that the melting 
point is: 
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a™ = -9.5. (56) 

The available 3D Monte Carlo simulations [26] unfortunately are not precise enough to 
provide an accurate melting point since the LLL scaling is violated and one gets values of 
a™ = —14.5, —13.2, —10.9 at magnetic fields 1, 2, 5T respectively. We found also that the 
theoretical magnetization calculated by using parameters given by ref. [26] is in a very good 
agreement with the Monte Carlo simulation result of ref. [26]. However the determination 
of melting temperature needs higher precision, and the sample size (~ 100 vortices) used in 
ref. [26] may be not large enough to give an accurate determination of the melting 
temperature (due to boundary effects, LLL scaling will be violated too). The situation in 
2D is better since the sample size is much larger. We performed similar calculation for the 
2D LLL GL liquid free energy, combined it with the earlier solid energy calculation [17,19] 

9soi = -7T7T + 2 log i-^ - — - 2.92. (57) 



2(3 A & 4tt 2 a, 



and find that the melting point a™ = —13.2. It is in good agreement with MC simulations 
[28]. 

Phenomenologically melting line can be located using Lindemann criterion or its more 
refined version using Debye - Waller factor. The more refined criterion is required since 
vortices are not point like. It was found numerically for Yukawa gas [60] that the Debye - 
Waller factor e~ 2W (ratio of the structure function at the second Bragg peak at melting to 
its value at T = 0) is about 60% at the melting point. Using methods of [62], one obtains 
for the 3D LLL GL model 



-2W 



0.59. (58) 



C. Fitting of the melting line. Values of the Ginzburg numbers of various 

superconductors 

In this subsection we use the above results to fit experimental melting line of several 
"fluctuating" superconductors. As an example on Fig. 2 of ref. [18] we presented the fitting 
of the melting line of fully oxidized Y Ba2Cu^0 7 [6]. Melting lines of two different samples 
of the optimally doped untwinned [2,38] near T c {Y Ba 2 Cuz0 7 -i), DyBa 2 Cuz0 7 [8] and 
(K, Ba)BiOs [39] are also fitted extremely well. The results of fitting are given in Table 2 
(To remind our convention, H c2 is defined as T c dH ^ T ^ \ t=t c rather than (often inaccessible) 
H c2 (T = 0)). 

Table 2 

Parameters of high T c superconductors deduced from the melting line. 



material 


T 

- 1 - C 


H c2 


Gi 


K 


7 


reference 


YBC0 7 _ s 


93.1 


167.5 


1.9 x 10~ 4 


48.5 


7.76 


[2] 




92.6 


190 


2 x 10" 4 


50 


8.3 


[38] 


YBC0 7 


88.2 


175.9 


7.0 x 10~ 5 


50 


4 


[6] 


DyBC0 6 . 7 


90.1 


163 


3.2 x 10~ 5 


33.77 


5.3 


[8] 


(K, Ba)Bi0 3 


31 


26 


5.3 x 10~ 5 


107 


1 


[39] 
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Our value for the Ginzburg number of YBCO and DyBCO estimated here are generally 
lower than the ones commonly believed in the literature. Often quoted value for YBCO is 
of order Gi = 0.01 (see page 1134 of commonly used ref.( [10]) ). The direct calculation 
from eq.(7) gives Gi = 0.003 for A = 1400 A, f = 15 A, and 7 = 7 (k = 93.3). Note 
however that these values are estimated from measurements at very low temperature. Our 
values of A and £ are fitted to the vortex physics experiments near T c and extrapolating 
using (admittedly questionable) two liquid model to T = give A = 931 A, £ = 18.7 A. 
Our values of dH f^ T ^ near T c are consistent with recent measurement [41] (about 2) and 
smaller than earlier ones. There is no consensus on values of k measured using the 
microwave technique at very low temperatures, however they are also generally smaller 
than 100 (smaller than 70 at T = and decreasing with temperature according to ref. [63] 
and from 50 to 60 in ref. [64]). This explains the difference of order of magnitude in Gi 
between the often used values and our fitting results (small k will lead a small Gi as 
Gi oc k 4 £ 2 T c 2 7 2 ). We emphasize that the actual small parameter in the theory is not Gi 
but rather 00 = \ / 2Gi7r 2 (see eq.(5)). Even for Ginzburg number as small as 2 x 10~ 4 this 
quantity is 0.2. As a result the effect of thermal fluctuations is important on a significant 
portion of the phase diagram. 

Recently it was found that thermal fluctuation are quite significant even in a low T c 
material (K, Ba)Bi0 3 . This is despite its lower critical temperature and very small 
anisotropy (and thereby very small Ginzburg number 5.3 x 10~ 5 ). Since this material is 
not a "strange metal" nor d - wave superconductor, its H c2 is directly accessible and there 
is no problem with direct estimate of Gi. ou — 0.1 for (K, Ba)BiOs is not much smaller 
than that of YBCO. There is therefore no surprise (contrary to a statement in ref.( [39]) 
that fluctuation effects are still experimentally observable in (K, Ba)Bi0 3 . In order to be 
able safely ignore thermal fluctuations the fluctuation parameter u should be of order 0.01 
in which case Gi should be smaller than 5 x 10~ 7 . These are the cases of most low T c 
materials. 



D. Magnetization jump at melting 

The scaled magnetization (of liquid or solid) is defined by: 



m{a T ) = --^—g(a T ) , (59) 
aaT 

while the LLL contribution to the magnetization is 

M LLL = ^- 2 ^m(a T ). (60) 

Using expressions eqs.(55) for solid and eqs. (51,53) for liquid the magnetization jump 
AM at the melting point a™ = —9.5 divided by the magnetization at the melting on the 
solid side is 

AM Am . . 

— - = = 0.018. 61 

M s m s 
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It is indeed small and is compared on Fig. 2 of ref. [18] (right inset) with experimental 
results of fully oxidized YBa 2 CusOj [6] (rhombs) and optimally doped untwinned 
Y Ba 2 Cu^O-j^s [4] (stars). The agreement is quite good. If the HLL contribution is 
significant (see next section) eq.(61) is expected to be violated. 



E. Specific heat jump at melting 

In addition to the delta function like spike at melting for specific heat experiments, the 
experiments also show specific heat jump. The theory allows us to quantitatively estimate 
it. 

The specific heat contribution due to the vortex matter is C = — Tj^-G^T, H). The 
normalized specific heat is defined as 

C 

c 



H 2 T 

where C m f = 4wk 2^ aT 2 is the mean field specific heat of solid. Substituting the definition of 
the scaled free energy eq.(15) and scaled temperature eq.(13a), we obtain: 

16/W bu \ 4/3 , 4/3 A / bu \ 2/3 , 

-^(2-2b + tfg"(a T ) 

Using our expressions for energy of liquid and solid we obtain the following specific heat 
jump at melting: 



Ac = 0.0075 ( 2 ^ + t ) 2 - 0.20G* 1 / 3 (b-l-t) (J^j 



2/3 

(62) 



Using the parameters of Y BCOj-s obtained by fitting the melting line, Table 2, we 
compare eq.(62) with the experimental result of ref. [2] in Fig. 2 of ref. [18] (right inset). 
Note that error bars are very large and also that disorder might be important [49], so that 
the agreement of the theoretical and experimental result of specific jump is not good as 
that of magnetization jump. 



V. HIGHER LANDAU LEVELS CONTRIBUTIONS. EFFECTIVE LLL MODEL. 

A. Where is the LLL approximation really valid? 

Contributions of HLL are important phenomenologically in two sections of the phase 
diagram. The first is at temperature above the mean field critical temperature T C (H) inside 
the liquid phase. The second is far below the melting point deep inside the solid phase. 
Naively in the solid phase, when "distance from the mean field transition line" is smaller 
than the "inter Landau level gap", 1 — t — b < 2b, one expects that higher Landau 
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harmonics can be neglected. More careful examination shows that a weaker condition 
1 — t — b < 12b should be used for a validity test of the LLL approximation [42] to 
calculate the mean field contributions in vortex solid. Additional factor 6 comes from the 
hexagonal symmetry of the lattice since contributions of higher Landau levels (HLL), first 
to fifth HLL do not appear in perturbative calculation. In the liquid state the question has 
been studied by Lawrie [65] using the Hartree - Fock (gaussian) approximation. The result 
was that the region of validity is limited, but quite wide, see Fig.l. 
In this section we will incorporate the leading HLL correction using Gaussian 
approximation and then compare the theoretical results with experimental magnetization 
curves. 



B. Gaussian Approximation in the liquid phase 



The free energy density beyond the LLL approximation is: 



G 



uoH, 



c2 



2llK 2 Vol 



log / Wexp / d 3 x-\d z i>\ 2 - a h \ij\ 2 + 



(63) 



where vol denotes volume. In the framework of the Gaussian (Hartree - Fock) 
approximation free energy is divided into an optimized quadratic part K, and a "small" 
part V. Then K is chosen in such a way that the gaussian energy is minimal. The gaussian 
energy is a rigorous lower bound on energy. Due to translational symmetry of the vortex 
liquid, an arbitrary U(l) symmetric quadratic part K has only one variational parameter 
e : 



K 



1 



J d 3 x^(\n^-b\^)+±\d^\ 2 + e\ 



The small perturbation is therefore: 

l~ = - / d 3 x 



(-a h - e) 



+ 2' 



The Gaussian energy consists of two parts. The first is the "Trace log" term: 



urn 



c2 



2ttk 2 voI 



log 



uHl b 

27TK 2 ^ 



Vnb 



(64) 



(65) 



(66) 



n=0 



The second is proportional to the expectation value of v in a solvable model defined by K 



2-KK 2 



2ttk 2 



E 



2y/2ir ^ Vnb + e 



+0J 



E 



2V2tt ^ Vnb + e 



(67) 



Both are divergent in the ultraviolet in a sense that at large n the sums diverge. 
Introducing an UV momentum cutoff which effectively limits the number of Landau levels 
to Nf — j — 1, the Trlog term diverges as: 
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I 00 I 

■= by \Jnb + e = ■= 

v n=0 v 



f A 3/2 + l £ _°_ A V 



(68) 



with the last term, function u being finite. The "bubble" integral diverges logarithmically: 



00 



1 



2V2tt ^ yjnb + e s/2n 



(69) 



where u' = -^-u(e,b). Substituting eq.(68) into the gaussian energy one obtains (in units 



1 2 



fi'Ga 



^3 



A 3/2 +a;( 1 A l/2 )2+ / 
V27T V 



•Oft 



2; v^TT 



A l/2 



(70) 



+2a;^^A 1/ V - eu + cj(w') 2 + «• 

V 27T 



The first term does not depend on parameters of the system and can be ignored (the 
renormalization of the reference energy density), while the second is u dependent and 
indicates that T c present inside ah is renormalized. Defining an = a r h + 2cj^ : A 1 / 2 , the 
above energy becomes: 



9Ga 



-UK 



fe Al/2 ' 2+ H"l)vb Al/2 



V2 

-eu' + uj(u') 2 + u. 



r I 

a h u 



(71) 



Thus the temperature T c and vacuum energy will be renormalized. The first two terms in 
free energy are divergent and linear in fluctuation temperature u, they will not contribute 
to any physical quantities like magnetization and specific heat. Minimizing the energy 
eq.(71), we get the gap equation: 



e = —a r h + 2oo u' 

Superscript "r" will be dropped later on. The function u(e, b) can be written in the 
following form 



(72) 



where 



v ( x ) = Yl 



n=0 



y/n + x - -(x + n + -)s + -(x + n 



_)2 

2 ! 



2. 1,3 

~3 {x -2 )2 - 



(73) 



(74) 



For the LLL model in the Gaussian approximation, v (x) = y/x. In the " Prange" limit [66] 
Gi — > 0, the free energy is 



2ttk 2 



V2ir V b J 



(75) 
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C. Integration of the HLL modes and the effective LLL model 



A method for treating HLL modes is integrating them and obtaining an effective LLL 
model. The (effective) LLL model is applicable in a surprisingly wide range of fields and 
temperatures determined by the condition that the relevant excitation energy e is much 
smaller than the gap between Landau levels b. Within the mean field approximation in the 
liquid e is a solution of the gap equation of eq.(72). For the LLL dominance region, we 
take a conservative condition e/e c = 1/20. One observes that, apart from the fields smaller 
than H LLL wO.l T for YBCO, the experimentally observed melting line and its 
neighborhood are well within the range of applicability of this approximation as shown in 
Fig.l. 

The effective LLL energy (we will use unit of energy density 2^in this subsection) 
functional is defined by: 

» oo 

9eff[H = log J J[V^V^e W {-f^ ,ro,A,^,]}, (76) 

where ipo is the LLL N = component field and the rest are denoted by ipi. Expanding the 
functional up to the fourth order in ip Q and to the second order in d z one obtains: 

g ef f[M = Ag + ^|^o| 2 + w IlllIM- (77) 



fLLlbPo] = ~ 
UJ 



2 

^\d z 1pHLL\ 2 ~ a h \i) H LL\ 2 + ^HLL\ 



The direct (no ip dependence) renormalization of energy is: 



/n 

J i=i 



Ag = ~ log / n^W;exp{-/^ L [^]}, (78) 



where the HLL energy is 

Ihll = — 



]j\d z i) H LL\ 2 - a h \tpHLL\ 2 + ^HLL\ 4 

where i/jhll = Y^i^Pi- To calculate Ag, we divide the fuLL into 



(79) 



Khll = I Q {W\ 2 - + \\0M 2 + ^iVf) (80) 

plus Jhll — K H ll- Taking e as the solution of the gap equation of eq.(72), Ag takes a 
form: 

Ag = geauss - 9 lll + 2 <|^ | 2 } «|^o| 2 } - <H 2 » ; (81) 
gLLL = log J Vtp Vtp*exp{-f LLL [ip }} 

Here geauss is the effective free energy of the full GL obtained in the first subsection of the 
current section, eq.(71), (\i>\ 2 ) is likewise the expectation value of \tp\ 2 in the full GL. The 
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quantity <?lll is the effective free energy calculated with variational parameter e and 
(l^ol 2 ) is the expectation value in the LLL GL. The consistency (or matching) requirement 
is: 



9eff = ~ 




This condition determines the value of At: 

At = 4 «H 2 ) - <|^ | 2 » = Auj [u'(e, b)] - 4 <|^ | 2 > (83) 




For YBCO, the correction At is small. The effective LLL GL approach achieves a 
simplification by starting from the LLL effective model with T c and other parameters 
renormalized to account for the contribution of the HLL modes. This is what we assumed 
in sections III and IV. In particular, this approach is very precise if we calculate the 
properties along the melting line. For example, the magnetization jump is mostly due to 
the fluctuation of the LLL modes, the background effective energy Ag will not contribute 
anything since it is the same on both sides of the melting line. 



D. The HLL contribution to Magnetization 

Generally when k is quite large and magnetization can be approximated by 

M = ~G(T,H). (84) 

The HLL correction will be calculated as follows. We numerically solve the gap equation 
(72) from which G{T,H) can be obtained. Then eq.(84) is used to calculate the 
magnetization of the full GL model in Gaussian approximation. The HLL correction is 
thus the magnetization of the full GL model in gaussian approximation minus the 
magnetization of the LLL contribution in gaussian approximation. We compare the 
experiments using following approximation. While the corrections due to HLL are 
calculated in gaussian approximation, the LLL contribution will be calculated 
nonperturbatively. The comparison of the theoretical predictions with the experiments for 
fully oxidized YBa2Cu^Oj [6], is shown on Fig. 3 of ref. [18]. We used the experimental 
asymmetry value 7 = 4 and values of T c , H C 2 and Gi from the fitting of the melting curve 
(see Table 2). The agreement is fair at intermediate magnetic fields, while at low magnetic 
fields is not good. It is expected that agreement is improved at higher fields. It is not clear 
whether magnetization (in contrast to magnetization jump at melting) will be strongly 
influenced by disorder, so at this time it is not possible to consider optimal doped YBCO 
magnetization curved more quantitatively. 

We comment that the theory of the full GL model (higher Landau levels included) beyond 
Gaussian approximation is required at low magnetic fields. Indeed experimentally it is 
often claimed that one can establish the LLL scaling for fields above 3 T for YBCO (see, 
for example, ref. [35]) as at low magnetic fields, the HLL contribution will be significant. 
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VI. SUMMARY 



The problem of calculating the fluctuations effects in the framework of the Ginzburg - 
Landau approach to vortex matter in type II superconductors is sufficiently precisely 
solved in the LLL approximation to allow quantitative description of the melting 
transition. We provided an evidence that metastable homogeneous state (the supercooled 
liquid state) exists down to zero fluctuation temperature by solving the large N Ginzburg - 
Landau model. Based on this understanding the supercooled liquid state is approached 
using methods of physics of critical phenomena (the Borel - Pade resummation technique). 
Applicability of the effective lowest Landau level model was subsequently discussed and 
corrections due to higher levels is calculated. 

The theory is then applied to quantitatively describe a great variety of experiments 
(confined to a region not far from T c ) including melting curves of 

YBCO, DyBCO, (K, Ba)Bi0 3 , magnetization curves, discontinuities of various quantities 
at melting. 

We speculate that any system of repelling objects (examples include classical one 
component plasma, electron gas...) exhibits similar features. The supercooled metastable 
state extends down to zero temperature. At this limit there is a well defined Madelung 
energy of the "ideal liquid". This ideal liquid is a pseudo critical point which controls the 
supercooled state possibly up to the melting temperature and might have universal 
features. 
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Figure captions 
Fig. 1 

Comparison of the experimental melting line for fully oxidized Y Ba2Cu^O-j [6] with our 
theoretical fitting. Applicabilty of the LLL approximation is between two lines, the solid 
LLL applicablity line and the (liquid) LLL dominace line. The GL model applicability line 
is also plotted. 

Fig. 2 

Free energy of solid (line) and liquid (dashed line) of the LargeN model as function of the 
fluctuation temperature l/|ar| 3 ^ 2 - The solid line ends at a point (dot) indicating the loss 
of metastability. 

Fig. 3 

Internal energy of the classical one component Coulomb plasma. The dashed line is the fit 
given by eq.(46) and the solid line is the fit given by eq.(44). 

Fig. 4 

The BP approximation for the free energy. BP3 and BP5 are the free energy results given 
by h 3 and h 5 . The dashed line Ti is the original perturbative expansion of order % in ref. 
[14] and the the dot dashed line % is the optimized expansion of order i. 
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Li and Rosenstein 
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